Mass Inflation in Brans-Dicke gravity 
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A detailed non-linear analysis of the internal structure of spherical, charged black holes that 
are accreting scalar matter is performed in the framework of the Brans-Dicke theory of gravity. 
We choose the lowest value of the Brans-Dicke parameter that is compatible with observational 
constraints. First, the homogeneous approximation is used. It indicates that mass inflation occurs 
and that the variations of the Brans-Dicke scalar inside the black hole, which could in principle be 
large in the absence of mass inflation, become small when mass inflation does occur. Then, a full non- 
linear numerical study of the black hole interior perturbed by a self-gravitating massless uncharged 
scalar-field is performed. We use an algorithm with adaptive mesh refinement capabilities. In this 
way, the changes in the internal structure of the black hole caused by mass inflation are determined, 
as well as the induced variations of the Brans-Dicke scalar, confirming, qualitatively, the indications 
given by the homogeneous approximation. 



I. INTRODUCTION 

It has long been known that the internal structure of 
an eternal Reissner-Nordstrom or Kerr black hole must 
be drastically altered by accretion. Firstly, a Reissner- 
Nordstrom or Kerr black hole has an inner horizon cover- 
ing a gravitationally repulsive singularity, which is qual- 
itatively different from a Schwarzschild black hole; it 
is unreasonable that the accretion of an infinitesimal 
amount of charged or spinning matter could lead to such 
a dramatic change in the black hole interior. Secondly, 
there are arguments, of both physical [l[ and mathemat- 
ical nature [2[, indicating that the Cauchy horizon (i.e. 
inner horizon) of the eternal charged or rotating hole is 
unstable against small (linear) perturbations. The natu- 
ral question is then, what is the endpoint of the instabil- 
ity? Answering it should shed some light on the internal 
structure of realistic astrophysical Kerr black holes. 

Studying general perturbations of a stationary axi- 
symmetric spacetime is a hard problem in general rela- 
tivity [3]. Thus, the aforementioned questions can be ini- 
tially considered in the context of spherical perturbations 
of the Reissner-Nordstrom black hole, as a simple model 
for the more challenging realistic case of the Kerr metric. 
This was the route followed by Poisson and Israel, in their 
seminal work where the phenomenon of mass inflation 
was unveiled 0, [H] : if ingoing and outgoing streams are 
simultaneously present near the inner horizon, then rel- 
ativistic counter-streaming between those streams leads 
to exponential growth of gauge-invariant measures such 
as the interior (Misner-Sharpe [(|) mass, the center-of- 
mass energy density, or curvature scalar invariants. Since 



'Electronic address: ppavelin@fc.up.pt 



^Electronic address: Andrew.Hamilton@colorado.edu 

fcTup.pt | 



* Electronic address: [crherdei 



this effect is causally disconnected from any external ob- 
servers, the mass of the black hole measured by an out- 
side observer remains unchanged by the mass inflation 
going on in the interior. But this inflation phenomenon 
causes the spacetime curvature to grow to Planckian val- 
ues on a spacelike hypersurface in the neighbourhood of 
the Cauchy horizon. 

A lot of effort has been put in understanding detailed 
aspects of the mass inflation phenomenon, within gen- 
eral relativity. In their original work, Poisson and Israel 
gave analytical arguments to support that the exponen- 
tial growth of the internal effective mass was sourced by 
the existence of both ingoing and outgoing fluxes - a rel- 
ativistic counter-streaming. According to these authors, 
the ingoing flux is infinitely blueshifted as it approaches 
the Cauchy horizon. But by itself this does not suffice 
to trigger mass inflation; an outgoing flux that induces a 
separation between the Cauchy and the inner apparent 
horizon is also necessary. 

Various other numerical and analytical studies in the 
framework of general relativity have shown that mass 
inflation develops naturally when a massless scalar field 
is accreted by a spherical charged black hole @, H, H, EH 
[Til . Il2l . [l3| . Perturbations of the scalar field propagate at 
the speed of light, and outgoing fluxes are produced as a 
result of the scattering of the ingoing scalar field by the 
space-time curvature, leading again to counter streaming 
and to mass inflation. Yet another example has been 
provided in [lil fl5j by studying the self-similar accretion 
of charged baryons into spherical black holes, showing 
that mass inflation occurs if an additional neutral matter 
component is present. In this case it is the relativistic 
counter-streaming between ingoing neutral matter and 
outgoing charged baryons that appears to be responsible 
for mass inflation. The physics of the relativistic counter- 
streaming that drives mass inflation inside black holes is 
further analysed in [l6j |. 

But it has not been explored, so far, if these dramatic 
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events taking place inside a black hole still occur in a 
more general theory of gravity. This is the programme 
we wish to start with the present paper. 

To be concrete, we shall consider one of the simplest 
extensions of general relativity, the Brans-Dicke (BD) 
theory of gravity [13], a well known example of a scalar- 
tensor theory in which the gravitational interaction is 
mediated by a scalar field, in addition to the metric. BD 
theory is parametrized by one additional parameter, u>, 
as compared to general relativity, which is obtained in 
the limit u> — ► oo. Hence, this theory passes all available 
experimental tests 18] if ui is large enough. The most 
stringent limits were obtained from precision measure- 
ments of the timing of signals from the Cassini mission 
which constrained u to be larger than 40000 [l^| . 

The key observation for the present study is quite sim- 
ple. If BD theory is minimally coupled to the Maxwell 
field, such a theory admits as a solution the standard 
Reissner-Nordstrom black hole of Einstein-Maxwell the- 
ory, irrespective of the value of u>. Indeed, the variations 
of the BD scalar field, <I>, are sourced by the trace of the 
matter energy-momentum tensor T, 



8ttT 
3 + 2w 



(1) 



which is zero for a classical electromagnetic field. Thus, 
the BD scalar $ can be taken as constant, and the theory 
reduces to Einstein-Maxwell with an extra (irrelevant) 
scalar field. 

Matter perturbations with non-vanishing energy- 
momentum trace will, however, induce variations of the 
BD scalar. Since this scalar determines the gravitational 
coupling to matter, it influences a back reaction of mat- 
ter on the geometry. A priori, there are reasons to ex- 
pect these variations to be large. The first reason is that 
equation |T]) with negligible or no sources, specialised to 
a stationary spherically symmetric spacetime in spherical 
coordinates, 



ds = gu(r)dt 2 + g rr (r)dr 2 + r 2 dVL2 



(2) 




yields, in a region where the metric coefficients have fast 
variations (this follows from (|54|) below) 



(3) 



thus, near a horizon, in particular near the Cauchy hori- 
zon, one could expect large variations of $ if an (even 
tiny) source term or perturbation in <E> exists. The sec- 
ond reason is that mass inflation could potentially lead 
to a large right hand side of ([1]). Thus, two natural (and 
related) questions follow: (a) does mass inflation still oc- 
cur? (b) are there significant variations of the BD scalar, 
and hence of the gravitational coupling to matter? 

In this paper we shall answer the above two questions, 
by studying the space-time geometry and the dynamics 
of the matter fields in the interior of an accreting spher- 
ical, charged black hole in the framework of BD gravity. 



We shall see that, choosing the lowest value of lu compat- 
ible with observational constraints, mass inflation occurs 
much in the same way as in general relativity, and no sub- 
stantial variations of the Brans-Dicke scalar are found. 
That $ varies little is actually a surprising conclusion, in 
view of the above arguments. As we shall argue, it is ac- 
tually the mass inflation phenomenon that, in our model, 
prevents larger $ variations from being produced. 

This paper is organised as follows. We start by in- 
troducing our model in section [TTJ In this section we 
discuss the equations of motion, assuming spherical sym- 
metry, specify the initial conditions for a full numerical 
integration and discuss the numerical algorithm used. In 
section |ni] we describe the homogeneous approximation, 
from which we extract both numerical and analytical 
information about the occurence of mass inflation and 
the variations of the Brans-Dicke scalar. In section IIVI 
we perform a full non-linear numerical evolution of the 
black hole interior perturbed by a self-gravitating mass- 
less scalar-field, using both a compact and a non-compact 
perturbation. In section [V] we draw conclusions. 



II. 



MODEL 



A. Field equations with spherical symmetry 



Consider the Brans-Dicke (BD) action: 



S = 



<J> <J>^ 



16tt£ 



M 



(4) 



where R is the Ricci scalar, $ is a real scalar field and u> 
is the Brans-Dicke parameter for which, throughout this 
paper, we shall take the minimal value consistent with 
observational constraints fl9j , namely 



uj = 40000 



(5) 



Lm is the matter Lagrangian, which we will assume to be 
the sum of the usual Maxwell contribution, Lp, and the 
contribution from a self-gravitating real massless scalar 
field tp: 

F 2 1 

£&r = £f + £ v = ~77 7T l P,a<P' a > (6) 

where F 2 — F a pF a P and F a p is the Maxwell tensor. 
Our signature choice is — h H — h The energy- momentum 
tensor of the matter fields is given by M T jiV — F T jJiV + 
T/xi/j where 



T 



-L II. 



J_ /„ „ „ 1 

47T 
1 

to\ r " tr "' 2 
The energy-momentum of the BD field is 



F F a 



1 



(7) 
(8) 



(9) 
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The equations of motion derived from ((4]) with ((6|) are 
the gravitational equations 



r ■- —( m t +*t ^ 



(10) 



equation Q for the BD field, with T = M T a a , and the 
matter equations 



□<p = , d*F = , 



(11) 



where * denotes Hodge dual. 

To solve the coupled system (TTJ, (fTU| and (fTT|) . we 
take a spherically symmetric ansatz written in double- 
null coordinates. Thus 



ds 2 



-2e 2a{u ' v) dudv + r 2 (u,v)dn 2 , (12) 



$ = $(tt,u) , F = F uv (u,v)duAdv , = <p(u,u) , (13) 



s 9 2 e 2CT 2 2 
- 2r(r,„„ + w,™) = + — r tp lU tp,v 

2 

+ J (»",,.$,„ + + y (w$,„$,„ + 2$$ „„) . (20) 

Equations (fT7|) - (fT5)) are two constraint equations, 
whereas equations (fTi?)) - (l2"2"]) are evolution equations. In 
the former equations, the contribution from the energy- 
momentum of the scalar field is 



(W) 2 

47T 



T v = 



4"7T 



(21) 



respectively, which represent, physically, the flux of the 
scalar field through surfaces of constant v and u, i.e. out- 
flux and influx. We shall impose that, initially, only in- 
flux exists, equation (|36[) . However, outflux is inevitably 
produced by scattering off the spacetime geometry. 

Finally, the wave equation (TT]) for the Brans-Dicke 
scalar field <I> is 



where u and v are taken to be ingoing and outgoing re- 
spectively. Under a gauge transformation of the form 
u — > function of u and v — > function of v, the coordi- 
nates it and v still preserve their null character; this will 
be of use later. 

The Maxwell equations (TTTj) are simply solved to yield 



F r 2 e- 2a 



constant = q , 



(14) 



r 



2<P,u<P,v 

3 + 2lu 



(22) 



B. Initial conditions 



Without loss of generality, we take the initial unper- 
turbed RN black hole to have mass 



where 



(15) 



is the electric charge. The electric field is therefore purely 
radial, as expected from spherical symmetry. The scalar 
field equation [[lip gives 



r 



(16) 



The uu, vv, uv and 99 components of the Einstein 
equations (|10j) . which are the only distinct non-trivial 
ones, give rise to the following equations (the left hand 
sides these equations are components of the Einstein ten- 
sor Gfu,)'. 



4r, u a u - 2r uu 2(^, u ) 2 + $ >uu - 2$, u a u 
' 



(17) 



$2 



4r v a v -2r vv 2(ip v y + $ vv -2$ v a v 
— — = — k 1" 



c + 2r„r„ + 2rr 



„2 _2ct (f, 



r 4 <J> <J> 

2(r,.^,»+r„$,») 
r <I> 



(18) 



(19) 



too = 1 



(23) 



which is sim ply a choice of units. Observational evidence 
[20l I21I [22I |23| suggests that astronomical black holes 
have large spins in several cases. We therefore take the 
unperturbed RN black hole to have a charge near the 
extremal value of unity, 



q = 0.95 . 

The corresponding radii r_ and r_| 
inner and outer horizons are 

r_ ~ 0.69 , r+ - 



(24) 

of the unperturbed 



1.31 



(25) 



Note that the mass of the black hole will change with 
accretion, but not its charge. 

The numerical integration is performed in double-null 
coordinates u, v over an integration box u = < u < 30, 
v = 5 < v < 20. Initial conditions are prescribed along 
the two null segments u — uq and v = vq. Prior to v = Vo, 
the geometry is Reissner-Nordstrom. We consider two 
cases, one in which the black hole is fed with a compact 
influx of scalar field <p that extends over a finite time, and 
another in which the black hole is fed with a non-compact 
influx that starts at a certain time and then continues 
indefinitely. Fig.[T]illustrates the initial conditions for the 
case of the compact influx; the case of the non-compact 
influx is similar. 
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V = Vq 




u — uo = 



FIG. 1: Initial conditions in the numerical integration box, 
defined by uo = < u < 30, vo — 5 < v < 20, for the case 
where the black hole is fed by a compact influx of scalar field 
ip. The unperturbed black hole has mass and charge mo = 1, 
q — 0.95. The influx along the null segment u — uo starts at 
vo = 5, has amplitude A = 0.05 or A = 0.1, and endures for 
Av = 1. The influx, depicted by a curve, is non-zero only for 
vo < v < vq + Av. The scale in the u and v axes is different to 
make closer contact with the Carter-Penrose diagrams shown 
in Figs. O and H 



The gauge freedom associated with the transforma- 
tions u — > f(u) and v — > f(v) can be used to make r 
linear in u and v on the initial null segments, so that 

r(u ,v) = v, r(u,v ) = v + ur jU (u ,v ) ■ (26) 

Defining r = r(uoj^o)i equations (|26p show that 

r = vo ■ (27) 

We choose 

ro = 5, (28) 

which is well outside the outer horizon. This ensures that 
the null segment u = u along which the initial influx ([55)) 
or (|4"Tj) of scalar field ip is specified is everywhere well out- 
side the horizon. The Misner-Sharpe [6J mass function, 
the total effective mass inside a sphere of radius r(u,v), 
is equal to 



m(u, v) 



1 



r 1 



q +4 r,„ 



2e 2 



(29) 



Equation (|2T)|) implies that r :U (u ,v ) in equation [[25]) is 
related to the initial mass mo and charge q of the black 
hole by 



r, u (uo,v ) 



1 


' 2 




- 1 
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(30) 



For the Brans-Dicke field $, which is observation- 
ally constrained to vary little outside the black hole, we 
choose initial conditions 



®{u ,v) = 1 = $(u,v ) . 



(31) 



This implies 



®,u(u, Vo) = = $,uu(W) V ) , 



(32) 



$,v(uo,v) = = 3>, vv (u ,v) . (33) 
The constraint equations (fTT|) and 1|18|) then become 



o-,u(u,w ) 



ro + ur, u (u ,v ) . , 



2r u (u ,v Q ) 



<r,v(uo,v) = -{tp, v (u ,v)} 2 



ip, u (u,v )r , (34) 
(35) 



We assume that there is no outflux along the initial 
null segment v — vo 



(36) 



Equation (|36[) implies that the scalar field (p is constant 
along the null segment v — Vq, and, since ip is denned 
only up to an overall constant, without loss of generality 
we take this constant to be zero, 



tp(u,v ) = 



(37) 



Equations p4)l and (f36|) together imply that <7 u (u, v n ) — 
0, which implies that a is constant along the null segment 
v = Vq. The constant can be chosen arbitrarily, since the 
choice (l2l)t still leaves a gauge freedom in the overall 
scaling of u, and we choose 



<t(u,Vq) = a-(u ,v ) = -ln(2)/2 , 



(38) 



equivalent to choosing g U v(uo,Vo) = 1 in the line- 
element (JT5J). 

The last initial condition to be specified is the in- 
flux tp tV (uo,v) of the scalar field along the null segment 
u = uo- We consider two cases, one consisting of a com- 
pact pulse of influx, the other a non-compact influx. The 
compact pulse represents an accretion event that takes 
place over a finite time, while the non-compact influx 
represents accretion that continues into the indefinite fu- 
ture. 

For the compact pulse, we choose the influx of scalar 
field along the initial null segment u = uo to be 

<p,v(u ,v) = A sin 2 i v <v<vq+Av, 

(39) 

vanishing outside the interval vo to vo + Av. We choose 
an interval 



Av = 1 . 



(40) 



The pulse has amplitude A. For A = 0, there is no 
pulse, and numerical integration yields the unperturbed 
RN black hole (see subsection III D|) . The behavior of u 
along the null segment u = uo follows from integrating 
equation (|33|) . 

For the non-compact influx, we choose the influx of 
scalar field along the initial null segment u = uq to be 



<p, v (uo,v) = — , V > v 
v 



(41) 



FIG. 2: The radial coordinate r as a function of u and v close 
to the inner ingoing apparent horizon (top) and outer ingo- 
ing apparent horizon (bottom), for the unperturbed RN black 
hole with m = 1 and q — 0.95. Horizons are at r_ ~ 0.6878 
and r+ « 1.31, in agreement with analytical expectation. 
Contours are plotted only for radii in the range of the color 
bar, r G [0.6876, 0.6879] in the top panel, r G [1.2, 1.4] in the 
bottom panel. 

C. Numerical algorithm 

We integrate the field equations using a second order 
algorithm similar to that used by other authors (see for 
example |13| ] and references therein). The four evolution 
equations can be written as 

YL = F S (Y k ,Y*,Y%) , k = l..A, (42) 

with Y J = r,a,Q,(p for j = 1 to 4 respectively. The nu- 
merical code determines the values of the variables Y^ in 
the (it, v) numerical grid solving, for each v = vq + n v Av, 
equation (|42|) for u — UQ + n u Au (starting at n u = 1) un- 
til the last point {v,f,vt) is finally reached. We performed 
standard numerical tests by running various simulations 
with different resolutions (Am and Aw), verifying that 
the code is indeed accurate to second order in Au and 
Aw. 

Although in double-null coordinates all the fields are 
regular, there is a numerical difficulty across the event 
horizon of the black hole, since r u diverges along the 



FIG. 3: Carter-Penrose diagram for the eternal RN black hole. 
Some constant r lines are depicted. The dotted line encloses 
a rectangular region in the u, v plane, where the integration 
takes place (in the chosen coordinates Uq = < u < 30, 
vo — 5 < v < 20). The dashed lines enclose two regions 
corresponding to the regions enclosed by dashed lines seen in 
Fig. [2] The same patterns of r = constant lines are seen. 

event horizon (see Fig. [5]- bottom). In order to maintain 
numerical precision, we used an algorithm with adaptive 
mesh refinement (AMR) in one direction. Given that 
strong gradients appear mainly in the it-direction, the 
AMR algorithm modifies the increment Au in such a way 
that Ar/r < 10~ 4 . This unidirectional approach is pos- 
sible since Courant-like stability constraints on the step 
size are avoided by the use of null coordinates. 



D. Unperturbed black hole 

Fundamental to the interior structure of a black hole 
is the location of its apparent horizons. Fig. [2] reveals 
the location of these horizons in the u, v plane for the 
unperturbed case. Fig. [5] shows the coordinate r as a 
function of u and v for the unperturbed RN black hole 
(without accretion, A = 0). Recall that an apparent 
horizon may be loosely defined as the boundary between 
the region where outwardly directed light rays (i.e. with 
u = constant and increasing v in null coordinates) move 
outwards (i.e. to larger r) and the region where outwardly 
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directed light rays move inwards. Apparent horizons are 
defined locally, and their locus is not, in general, a co- 
variant statement; it depends on the chosen time folia- 
tion of the spacetime. An illustrative example is that the 
Schwarzschild solution may be foliated such that it has 
no apparent horizon [241 ] , despite the presence of an event 
horizon, which is a covariant (and global) statement. 

In terms of the chart (|12p . apparent horizons are the 
set of points (u a ,v a ) such that 



dr 

T^(u a ,v a ) = . 



(43) 



This equation defines a curve in the u, v plane. Fig. [2] 
shows that the apparent horizons lie along u a — constant, 
and that r(u a ,v a ) ~ r± to good accuracy, meaning that 
the radii of the apparent horizons coincide with the ex- 
pected radii of the outer and inner event horizons, equa- 
tions ([25]). This is a good test of the code. As a further 
test we were also able to reproduce the results reported 
by E3. 

The content of Fig. [2] can also be understood in the 
context of the Carter-Penrose diagram of an eternal RN 
black hole. Fig. shows the Carter-Penrose diagram, 
outlining the regions covered by Fig. [21 The diagram 
reveals the same pattern of lines of constant r as Fig. [2] 



HOMOGENEOUS APPROXIMATION 



III. 



Results will be presented in section iTVl but first it is 
useful to gain i nsig ht by applying the homogeneous ap- 
proximation [9l. l25l l26j|. 

If the black hole is accreting at a low rate, it is rea- 
sonable to expect that the geometry, at least outside 
the region where mass inflation is occurring, would ap- 
proximate that of a stationary black hole, a Reissner- 
Nordstrom black hole. A stationary black hole has the 
property of time translation invariance 

|-0. (44) 

Inside the horizon, the time coordinate t is spacelike, and 
d, HH, H(| has called the proposition that equation l[l4")) 
holds inside the horizon the homogeneous approximation. 
fill point out that the homogeneous approximation (|44|) 
fails in the inflationary zone even in the limit of infinites- 
imal accretion rates. Rather, the homogeneous approxi- 
mation (|44p is equivalent to the assumption of symmet- 
rically equal ingoing and outgoing fluxes. The initial 
conditions prescribed in subsection IIIBI are asymmetric 
between ingoing and outgoing, so the homogeneous ap- 
proximation (|44| is not expected to hold. 

However, the homogeneous approximation (|44|) simpli- 
fies the mathematics, because it means that all quantities 
are functions only of the radial coordinate r. Although 
the assumption of symmetrically equal ingoing and out- 
going fluxes is not realistic, nevertheless the behaviour in 
this case may provide insight into the more realistic case 
of unequal ingoing and outgoing fluxes. 



A. Field equations 

The spherical symmetric, homogeneous line-element is 
given by eq. @. The tt and rr components of the Ein- 
stein equations (fTTH) are then 



2 1 

9rr 9rr ^9ri 

r £ g rr 



= ^ ("2f +* 2?) , (45) 
_(^T;+*T 1 T) , (46) 



9tt - 9tt9rr + rg' tt _ 8tt 



r 9tt9rr 

where the prime denotes a derivative with respect to 
r. Energy-momenta in angular directions follow from 
the energy-momenta in t and r directions, together 
with the equations of covariant conservation of energy- 
momentum. The non- vanishing components of T M „, 
F T A ,„ and *T M j, are, in {t, r, 8, 0} coordinates, 



IY2 



8?rr 4 



diag(l, -1,1,1) , 



diag(l, 1,-1,-1) 



1 




/ cj rr 


$' 


8ng rr 




\ Qrr 


$' 




9tt + 


$' 

LO 


167rg rr 


(- 


gu 





(47) 




(48) 


J)" 4 "" 


,(49) 


)' 


(50) 



whence 



M 



T 



!Y2 



4irg r 



(51) 



In the above we have taken into account that the black 
hole charge will give rise to an electrostatic field which is 
purely radial and consequently the only non-zero compo- 
nents of the electromagnetic tensor, F^ u are F tr = —F rt . 
It is then straightforward to show from Maxwell's equa- 
tions that 



Ftr — — 



qV-9tt.gr 



(52) 



The wave equation (jlip for the scalar field <p can be 
written as a first order equation 



9tt 



1 



(53) 



while the wave equation HJ for the Brans-Dicke field $ 
becomes 



£1 



9'tt 



9rr 



9tt g rr 



= 9r 



' 3 + 2uj 



(54) 



^From equations (|5Tj) and (|53[) - (p>4"l) , one reads that the 
source term for the variations of $ is proportional to 



9r 



1 



3 + 2u> g tt r 2 



(55) 



Thus, we conclude that the source term for $ variations, 
depends strongly on the ratio g rr / gtt- We shall return to 
this point in section Ull El 
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B. Generalized scalar field 

A massless scalar field has sound speed equal to the 
speed of light. As originally envisaged by [j], mass in- 
flation is expected to result from relativistic counter- 
streaming between an influx and outflux of scalar field. 
We now introduce a generalized scalar field, which allows 
us to vary the speed of sound and hence have models with 
and without mass inflation. Our goal is to see the conse- 
quences, within the homogeneous approximation, for the 
variations of the Brans-Dicke scalar. 

Consider a generalized scalar field ip whose action is 



,(<£>, X)d A x 



(56) 



where X = ip >fi cp ttl /2. It follows that the energy- 
momentum tensor of the generalized scalar field is that 
of a perfect fluid [27l ] 



•PT? = (j> + p)U»U v + P 8ji 



with 4-velocity 



V /= 2X 



(57) 



(58) 



and proper energy density and pressure p — 2Xp t x — P 
and p = £j v . If <p „ is timelike then the scalar field (p is 
formally equivalent to a perfect fluid. Generally, (p tU need 
not be timelike, but it is timelike in the homogeneous case 
(cf. subsection IIII Cp . 

In the general casep = p(ip, X), the pressure cannot be 
expressed only in terms of p, since for a fixed X and p at a 
given event it will always be possible to have different p t x 
by changing the value of ip at that event. On the other 
hand, if p = p(X) then both p and p will be functions of 
X alone, and consequently the relation between p and p 
is established once X is fixed. If p oc X^ +w ^^ with 
w = constant, then 



p = wp . 

The corresponding sound speed c s is 



(59) 



(60) 



For example, the Lagrangian density L v 



X yields 



the usual massless scalar field model with equation of 
state p — p, the model discussed in the previous sub- 
section IIII Al As another example, oc X 2 describes 
an ultra-relativistic fluid with equation of state p = p/3. 
As a final example, L v oc X n as n — > oo describes a 
pressureless fluid (dust). 



C. Homogeneous generalized scalar field 

The homogeneous approximation implies that tp is a 
function of r only. The only non- vanishing component of 



the 4-velocity U v oc U v oc <p iV , equation (|58]) . is therefore 
U r . Since g rr < inside the horizon, U r is 



U' 



(61) 



where the =F sign is — for tp' > and + for ip' < 0. The 
4-velocity U v is timelike, and therefore the behaviour of 
the generalized scalar field is identical to that of a per- 
fect fluid. The energy-momentum v Tjf of the generalized 
scalar field is 

vj* = f T 8 = v T <l> =p = wp = - w v T r t ( 62 ) 

and the trace of the matter energy-momentum tensor 
becomes 



M 



T = p(3w - 1) 



(63) 



The proper energy density p in terms of the metric co- 
efficients follows from energy-momentum conservation, 
T^ v '.„ = 0, which implies 



HK ^f\^-\ (64) 



where pi, gtu and r, are all integration constants; pi is 
the density at the surface r — r^, where gu is gtu- 



D. Initial conditions 



Equations (|6"3"|) and ([M]) define the right hand side (|5o| 
of the BD scalar equation of motion, and will be used in 
subsection IIII El In solving the latter equation, we shall 
work in the limit where pi is small, and impose initial con- 
ditions inside the unperturbed outer horizon, since the 
homogeneous approximation, which requires symmetri- 
cally equal influx and outflux, fails outside the outer 
horizon. Moreover, we take the initial conditions in equa- 
tion (ED to be 



= 0.95r 4 



l 

Qrri 



2m q 2 



(65) 



E. Results for the homogeneous approximation 

We now present results for the homogeneous approxi- 
mation applied to a generalized scalar field inside a spher- 
ical charged black hole in Brans-Dicke theory. We con- 
sider two cases, one with w = 0, the other with w = 1. 

In the case with w = 0, mass inflation does not oc- 
cur and large variations of the Brans-Dicke scalar $ are 
present. By contrast, in the w = 1 case, mass inflation 
occurs and $ variations are small. We will argue that it 
is actually mass inflation that prevents large <!> variations 
in the w = 1 case. 

Within the homogeneous approximation, the task is to 
obtain the behaviour of the metric coefficients in i.e. 
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r 

FIG. 4: The evolution of \g rr \, \gtt\ (top, solid and dashed 
lines respectively) and $ — 1 (bottom) with r, in the case 
w = with pi = 10 -4 . This behaviour of \g rr \, \gtt\ indicates 
the presence of a horizon around r ~ 0.69. The evolution of 
$ — 1 indicates that the BD scalar diverges as the horizon is 
approached. 



9m gtt, by solving (|16")) . and of the Brans-Dicke 

scalar field <&, by solving (|54[) . The matter scalar field ip 
is given by equation (|53|) . 

Fig. 0] shows the result for the numerical integration of 
9m 9tt and $, taking w — 0. The result is interpreted 
as follows. Equation ([29")) shows that the interior mass is 
large if |<? rr | is small, so a signature of mass inflation is 
that |<7 rr | is becoming exponentially small. This is not the 
case here. Rather, Fig. [4] (top) shows that |<7 rr | diverges 
and \gtt\ vanishes, around r ~ 0.69. This is characteris- 
tic not of mass inflation, but rather of matter hitting the 
inner horizon. Fig. [4] (bottom) indicates that the Brans- 
Dicke scalar is diverging as this horizon is approached. 
This is the expected behaviour of the BD scalar near 
a horizon, given a small variation of this scalar, as dis- 
cussed in the introduction, cf. equation ([3]) . It should be 
noted that this expectation relies on the homogeneous 
approximation, which breaks down at the outer horizon 
and outside. 

Fig. [5] shows the result for g rr , g u and <&, taking w = 
1. This corresponds to the standard case of a massless 




0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 



FIG. 5: The evolution of |ffrr[,|fftt| (top, solid and dashed 
lines respectively) and $ — 1 (bottom) in the w = 1 case for 
various choices of pi (from right to left in top figure and from 
dotted to solid line in bottom figure, pi = 10 -4 , 2 x 10 -4 , 4 x 
10 _4 ,8 x 10 -4 ). This behaviour of |<?rr|, \gtt\ indicates the 
existence of mass inflation. A noticeable jump in <E> occurs at 
the beginning of mass inflation, but it is followed by a gentle 
variation of 



scalar field, the same as considered in section [TH and can 
therefore be used to anticipate the results of section HV Al 
The behaviour is quite different from the w = case. 
Fig. [5] (top) shows that both g tt and g rr decrease sharply 
(and proportionally) around r ~ 0.69. This indicates 
the presence of mass inflation. A counter-intuitive fact 
shown in this plot is that for smaller values of the initial 
energy density pi, mass inflation is more abrupt. 

Fig. [5] (bottom) shows the behaviour of the Brans- 
Dicke scalar Around r ~ 0.69, where mass inflation 
starts, there is a clear jump in the variation of <E>, which 
is followed by a gentle decrease in $ to smaller radius. 
This indicates that the variations of are quite small in 
the mass inflation region. 

Having described the behaviour obtained from numer- 
ical integration we now try to understand it analytically. 
The case of w — 1, where mass inflation occurs, indicates 
that the variations of $ are small (less than one percent); 
thus $ w 1. It follows that | $ T M „| < | M T MJ/ |, and con- 
sequently may be neglected in comparison to M T MI/ 
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FIG. 6: The evolution of |gtt|/|<7rr| for w — (dark solid line) 
and for w = 1 (same cases as exhibited in Fig. [S]) 



in equations (J45J) and (J46J) . In this approximation, the 
problem reduces to that of general relativity coupled to 
matter and a scalar field, without back-reaction from the 
BD scalar. Moreover, Fig. [S] shows that the variations 
of g rr and g tt are quite fast in the mass inflation region. 
Thus, eqs. (145|) and (|46p may be approximated by 



Irr 



TCj rr 

9'tt 
rgttg rr 



-87T T\ — &7TW T£ , 



lir M T r 



These two equations imply that 



9rr gtt 



-grr OC (g tt ) 



(66) 

(67) 
(68) 



Fig. [6] shows the ratio |fltt|/|flrr| for the scenario with 
w = 0, where mass inflation does not occur, and the one 
with w = 1, where it does. For the latter, the ratio is 
approximately constant in the mass inflation region, in 
accordance with (p5|) , even though each metric coefficient 
individually varies rapidly, as shown in Fig. [5j 

Equation (|68|) . for w — 1, indicates that mass inflation 
occurs. It follows from (l64l). (IBTll and (f68l) that 



(69) 



9r 



1 



in the mass inflation region. Then 



g rr oc exp 



Ar 

r 



(70) 



where Ar scales roughly with ^fpl. Thus, the larger pi, 
the larger the scale Ar for variations of g rr . This agrees 
with the behaviour observed in Fig. [5] (top): mass infla- 
tion is more abrupt for smaller pi. 

Similarly, within this approximation scheme, we can 
obtain some analytical information for the variation of 



the Brans-Dicke scalar $ with r, using eq. (|54|) 
equation may be written as 



This 



In 



I gtt I 

\9rr 



( fftt )(l+-)/2 



m, (7i) 



where f(r) is some function 
portant for this argument, 
occur, as in the w — case, 
and the source term diverge, 
metric coefficients seen in Fi 
gent behaviour of as seen 
occurs, as in the w = 1 case 
equation (|TTj) to be written 



of r, whose form is unim- 
If mass inflation docs not 
both the coefficient of $' 
from the behaviours of the 
?. |U This sources a diver- 
in Fig. |U If mass inflation 
then equation (|68|) allows 



2 



In 



\9tt 



(l-w) 



G 



C(gt 



\(«>-l)/2 



/(r) 



(72) 

where C is a constant. Hence w = 1 is a special case 
where the coefficient of $' and the source term become 
bounded as gtt — > 0, and therefore the variations of $ 
remain small in the mass inflation region. This is con- 
sistent with the behaviour observed in Fig. [5] (bottom) . 
Thus, one concludes that, for w = 1, it is the existence 
of mass inflation that prevents large $ variations in the 
region around r ~ 0.69. 

To summarise, the homogeneous approximation indi- 
cates that for w = 1, mass inflation should occur and 
that the induced variations of the BD scalar (and there- 
fore of the coupling constant) within the mass inflation 
region are small. In the next section, we confirm these 
indications by a full non-linear numerical evolution of the 
model of section Ql] 



IV. RESULTS 

In this section we present results for a massless scalar 
field tp accreting on to a spherical charge black hole in 
Brans-Dicke theory. We consider two cases, addressed 
in turn in subsections IIV Al and IIV B[ one in which the 
initial influx of scalar field is turned on for only a fi- 
nite time, equation (|39[) . and another in which the initial 
influx starts at some time and then persists into the in- 
definite future, equation (|4"Tj) . 

The equations to be solved are the four equations given 
by (0 an d (fT§|) - ((2"2"]) . which constitute a complete set of 
equations. The remaining Einstein equations, eq. (fT7|) 
and l|18p. are automatically satisfied by covariant conser- 
vation of energy-momentum, and can be used to check 
the accuracy of the numerical integration. 



A. Full non-linear numerical evolution for a 
compact influx 

In this subsection, the black hole is fed with a compact 
pulse of massless scalar field, equation ([59"]) . 



10 




25- 
23i- 



1.4 
1.38 
1.36 
1.34 
1.32 
1.3 
1.28 
1.26 
11.24 
1 1.22 
'l.2 



16 18 20 




0.65 


30r 


0.648 


29 


0.646 


28 


0.644 


27 


0.642 




0.64 


3 26 


0.638 


25 


0.636 


24: 


0.634 


23: 


0.632 




0.63 





p 



1.4 

1.38 

1.36 

1.34 

1.32 

1.3 

1.28 

1.26 

1.24 

1.22 

1.2 



FIG. 7: The radial coordinate r as a function of the null co- 
ordinates u and v close to the inner apparent horizon when 
a single compact pulse with amplitude A = 0.05 (top) or 
A = 0.1 (bottom) and during Av = 1 is fed on to an ini- 
tially unperturbed RN black hole with mo = 1 and q = 0.95. 
The radius r is no longer a monotonic function of u as in the 
unperturbed case (Fig. [2]- top). The locus of the inner appar- 
ent horizon is a curve that interpolates between the maxima 
dr/dv = of the r = constant curves. For the pulse with 
larger amplitude, the inner apparent horizon changes more 
abruptly with v. 



Fig. [7] shows the value of the radial coordinate r as a 
function of u and v close to the inner apparent horizon 
when a single compact pulse with amplitude A = 0.05 or 
A = 0.1, and duration Av = 1, is fed on to an initially 
unperturbed Reissner-Nordstrom black hole with m = 1 
and q = 0.95. There are substantial changes compared 
to the Reissner-Nordstrom solution shown in Fig. [2] For 
example, r is no longer a monotonic function of u. In 
particular, the variation seen in Fig. [7] implies, accord- 
ing to the definition in section III C| that the locus of the 
inner apparent horizon in the u, v plane is not given by 
u a = constant; rather u a increases with v a . According to 
equation (|43|) . the locus of the inner apparent horizon is 
described by a curve that interpolates between the max- 
ima of the constant r curves in Fig. [7J The value of r 
on these maxima reveals that as time increases (u and 



FIG. 8: This plot shows a detail of the radius r in the vinicity 
of the outer horizon for the pulse with amplitude A — 0.05 
(top) or A — 0.1 (bottom). The behaviour is similar to the 
case without accretion A — (Fig. [2]- bottom), which can be 
attributed to the fact that the accreted mass is small. The 
value of it at large v (in the future) differs noticeably from 
the zero accretion case. 



v increase), the r coordinate of the inner apparent hori- 
zon decreases. This is to be expected. Indeed, the radii 
of the outer and inner horizons of a RN black hole vary 
with mass as 

dr± 
dm 



1± 



(73) 



Consequently, increasing the mass of a Reissner- 
Nordstrom black hole with a fixed q leads to a larger 
r + and a smaller r_. Of course, the geometry of an ac- 
creting black hole is not the same as the RN solution, but 
the observed decrease in the radius of the inner apparent 
horizon is consistent with equation (|73[) . 

By the same token, accretion is expected to cause the 
outer apparent horizon to become slightly larger com- 
pared to the unperturbed RN solution. This expectation 
is confirmed in Fig. [51 which shows a detail around the 
outer apparent horizon for the same simulations as Fig. [7] 
Fig. [5] shows two relevant features. Firstly, whereas at 
early v times the plot is quite similar to that of the un- 
perturbed case (Fig. [5]- bottom), at late v times there is a 



.. IAH 




FIG. 9: Heuristic Carter-Penrose diagram for the perturbed 
RN black hole. The initial perturbation is the curve depicted 
at the beginning of the integration box. The diagram is drawn 
to include the features observed in Fig. [7] and [5] namely that 
u a increases/decreases with v a for the inner apparent horizon 
(IAH)/outer apparent horizon (OUH). One expects the cur- 
vature to blow up at the Cauchy horizon (CH), and to become 
Planckian on a spacelike surface before the CH. 



slight but clear decrease in the u a coordinate of the outer 
apparent horizon (clearly seen in Fig. [5]- bottom, but al- 
ready noticeable in Fig. [5]- top). Secondly, r(u a ,v a ) 
is at late v times slightly, but again clearly, larger than 
the unperturbed value 1.31, as expected. The behaviour 
just described can be summarised in the heuristic Carter 
Penrose diagram of Fig. [5] 

Fig. [TUl shows log 10 (m) as a function of u and v, ex- 
hibiting the mass inflation profile close to the Cauchy 
horizon of the black hole. These figures show an expo- 
nential increase of the mass function as v increases. The 
increase is consistent with the theorems of [28|, [29[ , which 
show that, for a compact pulse as here, the mass tends 
to infinity and a mathematical singularity is attained at 
the Cauchy horizon, where v — > oo. The expectation of 
a curvature singularity is supported by Fig. 111! which 
shows the scalar curvature of the spacetime: 

R=^[l + 2e- 2 ° (d u rd v r + r 2 d u d v a + 2rd u d v r)] . 

(74) 

The same term e 2a d u rd v r that appears in |29p and is 
responsible for the mass inflation phenomenon also ap- 
pears in the Ricci scalar. Thus, not surprisingly, there is 
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FIG. 10: log 10 (m) as a function of u and v for A = 0.05 (top) 
and A = 0.1 (bottom), showing that mass inflation occurs 
as v increases, approaching the Cauchy horizon of the black 
hole. The mass function is quite sensitive to A. 




FIG. 11: Lines of constant Ricci scalar (actually log 10 \R\) for 
A = 0.05. 



12 




V 



FIG. 12: $ - 1 as a function of u and v for A = 0.05. The 

variations of $ are small slowly fading away after the initial FIG. 13: Non-compact influx: log 10 (m) as a function of u 
collapse of the scalar field pulse. and v for the perturbation (|4ljl . with A = 0.1. The plot 

shows that mass inflation occurs close to the Cauchy horizon 

of the black hole. 



a strong growth of the curvature that accompanies the 
growth of mass, which is manifest in Fig. 111! 

Although one expects that the curvature will diverge 
only on the Cauchy horizon (at v — ► +oo), which is a 
null surface, it is important to remark that the curvature 
should reach Planckian values at finite values of v, on a 
spacelike surface. Beyond this point one expects quan- 
tum gravity corrections to become dominant effects, and 
therefore one cannot trust the geometry given by any 
classical theory of gravity. 

Let us remark that the Cauchy horizon of the eternal 
RN black hole is composed of two lightlike sheets, as can 
be seen in Fig. [3l one corresponding to v = +oo and the 
other to u — constant in our coordinates. The initial 
conditions imposed create an asymmetry between these 
two branches of the Cauchy horizon. These conditions 
are appropriate to see in one branch the mass inflation 
phenomenon and in the other branch the evolution of the 
inner horizon under accretion. Thus, an outgoing (u — 
constant) observer sees mass inflation (as v increases, cf. 
Fig. rrU|) , but an ingoing (v = constant) observer sees 
approximately constant mass, for large u. 

Finally, Fig. [T2l shows $ - 1 as a function of u and v. 
It shows that the process of mass inflation does not seem 
to drive rapid variations of the Brans-Dicke scalar. The 
largest variations of <!> are associated with the initial col- 
lapse of the scalar field pulse onto the black hole. These 
fluctuations are transient, slowly fading away after the 
initial collapse of a self-gravitating massless scalar field 
pulse. The amplitude of the fluctuations appears to in- 
crease again, slowly, towards the Cauchy horizon. But 
this growth tends to stop and <fr to become a function 
of r only, cf. Fig. [7] - top. This is consistent with the 
homogeneous approximation. 
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FIG. 14: Non-compact influx: log 10 \R\ as a function of u and 
v for the perturbation (|41[) . with A = 0.1. The behaviour 
of the curvature is qualitatively different from the compact 
pulse case (cf. Fig. I12[) . 



B. Full non-linear numerical evolution for a 
non-compact influx 

The model of the previous subsection, in which the 
black hole is fed with a compact pulse of scalar field over 
only a finite time, is not realistic. The model is equivalent 
to assuming that the black hole remains isolated forever 
after accreting the pulse of scalar field. Real black holes 
are however never completely isolated, and they will con- 
tinue to accrete at some rate or other until something 
completely different happens, like the black hole Hawk- 
ing evaporating, or the Universe coming to an end in a 
big crunch. 

In this subsection, the black hole is fed with an influx 
of massless scalar field that starts at some time but then 
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FIG. 15: Non-compact influx: $ — 1 for the perturbation (|41|) . 
with A = 0.1. The Brans-Dicke scalar variations are qualita- 
tively different from the compact pulse case (cf. Fig. Ill|l . 

continues into the indefinite future, equation (|4ip . 

The results of the simulation for the structure of the 
inner and outer horizon are qualitatively similar to those 
of the compact pulse. The initial conditions were chosen 
so as to have a similar black hole mass variation, of the 
order of 1% during the time spanned by the simulation, 
for both the compact and non-compact influx. In the 
latter case, accretion would continue to increase the black 
hole's mass. 

Figs. [13] to [15] show the mass function log 10 to, the 
Ricci scalar log 10 R, and the variation $ — 1 of the Brans- 
Dicke scalar. The behaviour of the mass function, in the 
integration domain, is similar to the compact pulse case. 
There is a subtle difference in that ingoing (v = constant) 
observers see more mass inflation than in the compact 
case. The curvature and Brans-Dicke scalar variations 
exhibit, on the other hand, a noticeable difference. For 
the non-compact influx they are monotonic functions of 
v, for large u, whereas for the compact pulse, they stop 
growing at some v and resume growing for larger v. The 
transition seems to be associated to the end of accretion. 



V. CONCLUSIONS 

We have studied the internal structure of spherical 
charged black holes in the framework of the Brans-Dicke 
theory of gravity, together with the variations of the 
Brans-Dicke scalar, which determines the gravitational 
coupling. 

The initial expectations drawn from naive arguments 
given in the introduction suggested possible large varia- 
tions of the Brans-Dicke scalar $ near the inner horizon of 
the black hole. The use of the homogeneous approxima- 
tion indicated that these expectations were correct: in a 



model with w = 0, the Brans-Dicke scalar had a divergent 
behaviour near the inner horizon. However, in a realistic 
scenario, accretion will generate counter-streaming and 
hence mass inflation. In turn, mass inflation changes the 
dynamics of the Brans-Dicke scalar and prevents large 
variations of this scalar from occurring. This conclusion 
was obtained by comparing the behaviour of models with 
w = and w = 1. We emphasise that the non-existence 
of mass inflation in a single-fluid model with w = is 
unrealistic, because even the tiniest admixture of a sec- 
ond fluid would allow rclativistic counter-streaming and 
therefore mass inflation to occur. 

One counter-intuitive result drawn from the homoge- 
neous approximation is that mass inflation is more abrupt 
for a smaller value of the initial perturbation, in agree- 
ment with the conclusions of (l6j . 

We have performed full non-linear numerical simula- 
tions of the evolution of the black hole perturbed by both 
a compact and a non-compact influx. In the case of the 
compact pulse, one observes some (expected) changes in 
the black hole interior, namely the growth of the outer 
horizon, decrease of one of the light- like sheets of the in- 
ner horizon and exponential growth of the effective mass 
function and curvature near the other light-like sheet of 
the inner horizon. In this case, the black hole settles 
down after accretion to a new stationary state, as illus- 
trated in the Carter-Penrose diagram in Fig. [5] The rel- 
ative variations of the Brans-Dicke scalar are small (be- 
tween 10 -6 and 10 -7 ) in the integration region when the 
compact influx increases the black hole mass by around 
1%. These variations grow slowly towards the Cauchy 
horizon, but seem to stabilize. The end of accretion is 
noticeable, in the case of a compact influx, in both the 
curvature and Brans-Dicke scalar variations. 

For a non-compact influx, on the other hand, the black 
hole does not settle down to a stationary state. This is 
a more realistic scenario, since real black holes are not 
expected to stay isolated forever, after an initial accretion 
period. For this type of perturbation, the changes in the 
black hole structure are similar to those of a compact 
influx. Noticeable qualitative difference appear, however, 
in the curvature and Brans-Dicke scalar variations. For 
an accreted mass of around 1% of the black hole mass (in 
the time spanned by the simulation) these variations are 
still small, albeit around one order of magnitude larger 
than for the compact influx. 
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